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We introduce the concept of nested topological order in a class of exact quantum lattice Hamilto- 
nian models with non-abelian discrete gauge symmetry. The topological order present in the models 
can be partially destroyed by introducing a gauge symmetry reduction mechanism. When symme- 
try is reduced in several islands only, this imposes boundary conditions to the rest of the system 
giving rise to topological ground state degeneracy. This degeneracy is related to the existence of 
topological fluxes in between islands or, alternatively, hidden charges at islands. Additionally, island 
deformations give rise to an extension of topological quantum computation beyond quasiparticles. 
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The concept of topological orders [l[ offers the possibil- 
ity of finding new states of matter with a common picture 
of string-net condensation Q and other variants thereof 
[H . They correspond to examples of long range entangle- 
ment in quantum many-body systems where those cor- 
relations emerge in quantum states that are encoded in 
non-local degrees of freedom of topologically ordered sys- 
tems. Their global properties are the source for yet an- 
other application as the suitable systems to implement 
topological quantum computation |1, H H 0, a form 
of fault-tolerant quantum computation intrinsically resis- 
tant to the debilitating effects of local noise. Quantum 
field theories with an spontaneous symmetry breaking 
mechanism of a continuous gauge group down to a dis- 
crete group have been pro posed as a scenario for realizing 
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their physics 

In this paper we introduce the concept of nested topo- 
logical order in a class of quantum lattice Hamiltonians. 
Our starting point are the family of Kitaev's models [1], 
which are labeled by a discrete gauge group. Such mod- 
els can be modified [J3| introducing an explicit symme- 
try breaking mechanism. Our aim is to study the ef- 
fect of 'nesting' subsystems with a reduced symmetry 
inside systems with the complete gauge symmetry. We 
will consider a topologically ordered system divided in 
two regions, say A and C, and show that it is possible 
to partially destroy the topological order in region C in 
such a way that this imposes boundary conditions to the 
subsystem A. The system C can take the form of several 
islands, which is why we talk about 'nested' topologi- 
cal order. The boundary conditions induce a topological 
ground state degeneracy which is due to the possible val- 
ues of certain fluxes in between islands. As we will see, 
the values of these fluxes correspond to the types of do- 
main walls that exist in C. If we allow the region C 
to be deformed, then islands can be initialized, braided 
and fused, giving an interesting extension of the ideas of 
topological quantum computation beyond quasiparticles. 

The models that we consider are string-net conden- 
sates in a 2D lattice 0. The configurations of the 
lattice are regarded as string-net states: a collection of la- 
beled strings meeting at branching points. A string-net 
is closed if certain conditions hold at branching points 



and there are no loose ends. The ground state is a su- 
perposition of all possible deformations of such closed 
string-nets, and excited states correspond to configura- 
tions with loose ends: quasiparticle excitations appear at 
the ends of strings. Now, to such system Hamiltonians 
we can add string tension terms, which penalize with a 
higher energy those configurations with longer strings. 
As such terms get more important with respect to the 
original ones, longer strings become less relevant in the 
ground state and finally the topological order is destroyed 
as excitations get confined. Alternatively, we can add 
suitable terms so that only part of the topological order 
is destroyed. This is in fact the case for the Hamiltoni- 
ans Hq' M that we consider (JTJ) , which are labeled with a 
discrete group G and two subgroups N C M C G, with 
N abelian and normal in G. If N = 1 and M = G, we 
have the original topologically ordered models with gauge 
group G considered by Kitaev Otherwise, the gauge 
symmetry is reduced the quotient group G' = M/N. In 
particular, if N = M the topological order is completely 
destroyed. 

Topological phases. The systems of interest are 
constructed from a two-dimensional oricntable lattice, of 
arbitrary shape. At every edge of the lattice we place 
a qudit, a |G|-dimcnsional quantum system with Hilbcrt 
space TL' G and a basis \g) labeled with the elements of G. 
The Hamiltonians read as follows 11511 
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where the sums run over the set of vertices V, faces F 
and edges E. Explicit expressions for the terms in ([I]) 
will be given below, but before that, we will discuss their 
physical content. First, all the terms are projectors and 
commute with each other, so that the ground state is 
described by conditions of the form P|GS) = |GS) with 
P either a vertex, face or edge operator. Excitations are 
gapped and localized; they correspond to violations of 
the previous conditions and so can be related to vertices, 
faces and edges; they are regarded respectively as electric, 
magnetic and domain wall excitations. 
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FIG. 1: Examples of lattice constructions. Although all the 
edges must be oriented, only the orientation of some of them 
is shown. The Ti,i = 1,2,3,4 are triangles; the light thick 
arrow shows their orientation, n and T4 are dual, the others 
are direct, a is a closed ribbon; the projectors K^' c give the 
charge in the region S that a encloses, p is an open ribbon; 
the projectors J&' T give the domain wall flux in the region 
T in the direction of the arrows, a and j3 are minimal closed 
ribbons, enclosing respectively a single vertex and face. 



related and the excitation types, labeled as (R,C), are 
dyons: C, the magnetic part, is a conjugacy class of G 
and R, the electric part, is an irrep of Np, the group 
N<7 := { g G G\grc = r c9}> where tq is some chosen 
element of C. These charges have a topological nature: 
if there are several excited spots in the system, far apart 
from each other, there exist certain global degrees of free- 
dom which cannot be accessed through local operators. 

In the general case Hq' M there are two new phenom- 
ena, quasiparticle condensation and the appearance of 
domain wall excitations. The latter have an energy pro- 
portional to their length and can be labeled by pairs 
(R,T), with T G M\G/M and R an induced repre- 
sentation in M of an irrep of the group :— { m G 
M I mr^Af = rr-M}, where is some chosen element 
of T. Thus there exists a flux related to domain walls, 
with values (i?, T); it is conserved in the absence of quasi- 
particle excitation, so that domain walls only can end at 
them. As for condensation, we will comment upon it 
below. 

Ribbon operators. In order to motivate the intro- 
duction of ribbon operators, we first note that dyons, the 
excitations of our system, are located at vertex-face pairs, 
which are called sites. In Fig. [1] sites are represented as 
dotted lines connecting the vertex to the center of the 
face. The basic connectors between sites are triangles: 
just as an edge connects two vertices, triangles connects 
two sites. A direct (dual) triangle r is composed by two 
sites and a direct (dual) edge e T , see Fig. [T] Triangles 
can be concatenated to form ribbons connecting distant 
sites. Ribbons are open if they connect disjoint sites and 
closed if their ends coincide. The point is that it is pos- 
sible to attach to each ribbon p certain operators F' p 1 ' 9 , 
h, g G G, which are very well suited to represent excited 
states. For example, any state with only two dyons is a 



linear combination of the states -F^IGS), with p any rib- 
bon connecting the sites where the dyons are located [1]. 
In fact, one can consider that ribbon operators represent 
a process in which a particle-antiparticle is created in one 
end of the ribbon and one of them is moved to the other 
end. 

In order to describe ribbon operators, we start with 
triangles, which arc the smallest ribbons. Recall that a 
triangle is formed by two sites and one edge, direct or 
dual. Triangle operators act on the qudit attached to 
that edge, and the action depends on the orientation of 
the edge and the type of the triangle. The four possi- 
ble cases are illustrated in Fig. Q] With the notation 
of that figure, we have = 5 gA J2 k \ hk )( k l F k § = 

\ 9 - 1 ){g- l \, F%> = \ g){g \ an d = <*„,! EfclM-m 
where the sums run over G. Then if p is a ribbon formed 
by the concatenation of the ribbons p± and P2, we set 
F^ a = y\ F^ k F^ lhk ' k ~ lg . The terms in the Hamilto- 
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nians (JT]) are built from ribbon operators. Let F^ v := 

\U\- 1 Y.»&jY, v &r F ? v for an y subgroups U,V C G. 
Then Af := F^ G , Bf := F} N , T e M := F T 1M and 
:= F^, G , with a and (3 suitable minimal closed rib- 
bons as in Fig. [T]and r (t') a direct (dual) triangle with 
e = e T . 

Ribbon operators commute with all the vertex oper- 
ators A G and face operators Bj, except with those at 
their ends. Moreover, they can be characterized by this 
property [15|. This suggests considering, for closed rib- 
bons a, those ribbon operators which commute with all 
vertex and face operators, so that they 'forget' the single 
end of fj. It turns out that a linear basis for such op- 
erators is given by a family of projectors K^' G , labeled 
with the charge types (R, C) of the system Hq. In fact, 
if a is a boundary ribbon, that is, a closed ribbon en- 
closing certain region S as in Fig. [lj then K^ G projects 
out those states with total topological charge (R, C) in 
S. As a result, the ground state of Hq can be described 
by the conditions 
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which must hold for all boundary ribbon a. This amounts 
to impose that all disc shaped regions must have trivial 
charge because K^ 1 = F G ^, where e is the identity repre- 
sentation. In systems with Hamiltonian Hq M we can use 
the projectors K^ G to describe condensation. Namely, 
for some charges |15| we have a ground state expectation 
value (K^ G ) > for any boundary ribbon a. showing 
that there exist a non-zero probability of finding such 
charges in a given region. 

Domain wall types can be obtained in a similar fash- 
ion in systems with Hamiltonian Hq M . For any open 
ribbons p, those ribbon operators that commute with all 
vertex operators A 1 ^ and face operators Bj are linear 
combinations of certain projectors J^' T , with (R,T) a 
domain wall type. If p crosses an area with domain wall 
excitations then Jjf' T projects out those states with total 
domain wall flux (R, T) across p. For example, in Fig. [T] 



p will measure the flux of the excited region T in the 
direction of the white arrows. 

The ground states of JT]) can also be described in terms 
of conditions for ribbon operators, in particular by 

F™ N \GS) = |GS), M |GS) = |GS), (3) 

where a and p are arbitrary boundary and open ribbons, 
respectively. The first condition is related to vertex and 
face excitations, and the second to edge excitations. 

Nested phases. We are now in position to discuss a 
more complicated system. In particular, we want to con- 
sider a surface divided in two regions of arbitrary shape, 
A and C, plus a third region B which is just a thick 
boundary separating them, included so that the Hamil- 
tonian does not have to change abruptly from A to C. 
The idea is to have a local Hamiltonian such that condi- 
tions ([21 are satisfied in A, conditions (J3j) in C and the 
conditions 

F? N \GS) = |GS), (4) 

with a an arbitrary boundary ribbon, in the whole sys- 
tem. The last condition is needed to ensure that domain 
wall flux is preserved through region B, a key ingredient 
of our construction as we will see. The ground state of the 
Hamiltonian Ho ■= — J2v ~ J2f Bf 1S described pre- 
cisely by (TJ|. In addition, Ho commutes with Hq, Hq . 
Indeed, a Hamiltonian of the form H' = Hq+XHo, A > 0, 
only differs from Hq in the gap for some excitations, and 
the same is true for Hq M . The Hamiltonian that we 
want to consider takes the form H = Ho + XHq + pHq M , 
where A, p > vary spatially so that A = 1 and p = 
in A and A = and p = 1 in C. If we take Xp = 0, the 
ground state has the desired properties but there exists 
some local degeneracy at B. This local degeneracy can 
be lifted if A and p are allowed to overlap, but on the 
other hand if the overlap is too big, it could produce a 
level crossing taking the ground state of H out of that of 
Ho, which spoils conditions 0. 

Quasiparticle dilution. Our aim is to understand the 
effects of the nested region C on the topologically ordered 
region A. A first effect is the possibility to locally create 
or destroy single quasiparticle excitations in the vicinity 
of the A-C border, something prohibited in systems with 
Hamiltonian Hq due to charge conservation. In terms 
of ribbon operators, this is reflected in the fact that for 
any p\ connecting C to A, as the one in Fig. [2ja), a 
state of the form £ meM F™< m f|GS), n e N, contains 
no excitation at C. In terms of quasiparticle processes, 
this corresponds to create a particle-antiparticlc pair in 
A and then move one of them into C , where it disappears 
because it is condensed. 

Domain wall dilution. A second effect is related to the 
existence of domain walls in region C. Consider again a 
ribbon p 2 connecting C to A, see Fig. HJa). Some of the 
states of the form = c h, g F h ' 9 \GS) , c Kg e C, 

will contain edge excitations all along the portion of p 2 
contained in C, for example those with Ch. g ^ for some 
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FIG. 2: In this figure regions A, B and C are shaded 
respectively with medium, dark and light gray. Ribbons 
Pi, i — 1, ... ,8 are displayed as pairs of solid and dashed par- 
allel lines which correspond respectively to their direct and 
dual edges. Light spots at the end of ribbons represent exci- 
tations in A and the dark one an excitation in C. The striped 
areas are domain wall excitations, (a) Due to condensation, 
suitable ribbon operators attached to pi will create an excita- 
tion in A but no excitation in C. Ribbon operators attached 
to p2 can create a domain wall excitation in C. The resulting 
state ip is such that J^ T \ip) = \tp). (b) Both ps and pa 
measure the flux in between the islands. If O is an operator 
with support in the shaded area an takes ground states to 
ground states, it cannot change the flux, (c) If the previous 
islands are deformed till they fuse, the flux measured by p7, 
ps will remain the same as it was for ps, pe- If it is nontriv- 
ial, opposite border charges are present at the sides of the 
meeting point. 



g S G,h $ M. These excitations form a domain wall, to 
which we can relate a type or flux given by the projector 
J^ c , where p^ is a ribbon that lies in C and crosses 
the domain wall, see Fig. HJa). Such a ribbon can be 
deformed without crossing any quasiparticle excitation 
onto another ribbon p^ that only has its endpoints in 
C and thus avoids the domain wall, so that Jp 3 ' T \^p) = 
J^' T \ip) due to ((3]). Both ribbon operators are measuring 
the same domain wall flux. However, in the case of p4 the 
flux is being measured in A, where the domain wall gets 
diluted as it turns into a condensed string. Note that 
Jrc cann °t detect changes in the interior of G. In this 
regard, if we restrict our attention to region A, domain 
wall flux projectors from ribbons like p4, that is, which 
enclose a portion of the A — C border, can be related to 
charges (R, T) that lie in that piece of the A — C border. 

Induced topological fluxes. Things get even more in- 
teresting if we consider that C consists of several disjoint 
parts. For example, consider a plane and choose as the 
region C two islands C\ and C 2 , see Fig. [Hb). Now in- 
stead of considering a domain wall flux coming out from 
a region of B (such as the one measured by p4 in Fig. 
[2{a)), we consider the flux in between the two islands 
(as indicated by the arrows in Fig. Hlb)). This is the 
flux measured by the projectors J^' T , where p 5 is any 
ribbon that connects the islands, as in Fig. (Hb). The 
point is that such a flux is a global (topological) prop- 
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erty as long as the islands arc distant. Indeed, measuring 
the flux requires an operator with a support connecting 
C\ and C2. And, if an operator changes the flux, its 
support must loop around C\ (or C2). Suppose to the 
contrary that O is an operator that leaves the ground 
state invariant and has a support not enclosing C\, as 
the shaded region in Fig. (Hb). Let p§ be another rib- 
bon connecting the islands but lying outside the support 
of O. Due to dl we have J* 5 < T |GS> = J*' T |GS>, so 

that [J£ T , 0]|GS) = [J^ T , 0]|GS) = and thus O does 
not change the flux. Those operators which do change 
the flux arc related to processes in which a particle- 
antiparticlc pair is created, one of them loops around C\ 
and they meet again to fuse into a charge that disappears 
into C\. 

Topologically protected subsystems. It follows that 
there exist a topological degeneracy in the ground state, 
related to the distinct values that the flux in between C\ 
and C2 can take. For example, if N = M = 1 the flux 
can take any value g G G. In general, for a C composed 
of multiple disconnected regions, the degeneracy of the 
ground state depends on N, M and the topology of A. 
Now, it is natural to ask how does this protected space 
compares with the one due to to the existence of several 
separated excitations in A. In other words, do islands add 
something new? This can be positively answered through 
an example: two excitations give no protected subspace 
[J] , but we have just seen the contrary for the case of two 
islands. Perhaps more dramatically, for abelian groups 
G the protected subsystem is always trivial whatever the 
amount of excitations, but this is not the case for islands. 
Nevertheless, islands can be compared to excitations, in 
the following sense. An island can hold certain charge 
values, which can be measured using ribbon operators 
that enclose the island, as in the case of an excitation. 
The difference between a charged island and a charged 
excitation is that the local degrees of freedom of the exci- 
tation become global in the case of the island: this is the 
origin of the additional dimensionality of the protected 
subsystem. 

Braiding. The physics of the system so far has a static 
nature. If we want to consider the setting as an sce- 



nario for quantum computation, then the possibility of 
dynamically deforming the region C must be included 
in it. Such deformations need not be strictly adiabatic, 
but the state should be kept in the subspace defined by 
conditions ([2]|4]) at all time. We can then braid islands 
to perform unitary operations, in complete analogy with 
quasiparticlc braiding. It is also natural to enrich the 
physics by considering islands with different (TV, M) la- 
bels, increasing the variety of protected subsystems. 
Fusion. We must consider also the analogue of the quasi- 
particle fusion processes, which is the way in which mea- 
surements arc carried out in topological quantum com- 
putation. There arc two natural ways in which global 
degrees of freedom can be made local. The first is to 
decrease the size of an island till it dissapears leaving 
a small charged region. The outcome of such a process 
is the charge, which can be measured but not changed 
locally. The second way is closer to the idea of fusion. 
Indeed, it is also a fusion, but of islands instead of quasi- 
particles. The idea is depicted in Fig. Etc). As two is- 
lands of the same (N, M) type get closer, some of the 
ribbon operators connecting them become small and thus 
the flux between the islands is exposed to local measure- 
ments. If we continue the approach till the islands meet, 
the flux will take the form of a domain wall excitation at 
the meeting place, as in Fig. [He). Due to confinement 
the domain wall can decay to several smaller walls, but 
there is something that will not disappear, the two bor- 
der charges in its ends on B. As explained in the caption 
of Fig. (He), the appearance of this border charges can 
be seen directly in terms of ribbon operators. Regarding 
the initialization of the system, reverse processes can be 
used. That is, if an island is divided in two, the topolog- 
ical flux in between them will be trivial, and if an island 
is created from the vacuum, it will have trivial charge. 
In both cases the reason is that topological properties 
cannot be changed by local processes. 
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